Dielectric quantification of conductivity limitations due to nanofiller size in 
conductive powders and nanocomposites 
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Conducting submicron particles are well-suited as filler particles in non-conducting polymer ma- 
trices to obtain a conducting composite with a low percolation threshold. Going to nanometer-sized 
filler particles imposes a restriction to the conductivity of the composite, due to the reduction of 
the density of states involved in the hopping process between the particles, compared to its value 
within the crystallites. We show how those microscopic parameters that govern the charge-transport 
processes across many decades of length scales, can accurately and consistently be determined by 
a range of dielectric-spectroscopy techniques from a few Hz to infrared frequencies. The method, 
which is suited for a variety of systems with restricted geometries, is applied to densely packed 7- 
nm-sized tin-oxide crystalline particles with various degree of antimony doping and the quantitative 
results unambiguously show the role of the nanocrystal charging energy in limiting the hopping 
process. 

PACS numbers: 73.22.-f, 72.80.Tm, 73.63.Bd, 77.84.Lf 
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I. INTRODUCTION 

Small submicron particles are well-suited as fillers in 
non-conducting polymer matrices to obtain a conducting 
composite with a percolation threshold (far) below 1 %. 
The low percolation threshold is due to the formation 
of airy aggregates of conducting particles, in which the 
particles are grown together by diffusion-limited cluster 
aggregation, creating a network with a fractal dimension 
around 1.7 These airy aggregates can be thought 
of as conducting spheres forming a 3-dimensional per- 
colating network around the expected aggregate filling 
fraction 0.16. As a consequence of the fractal structure 
within the aggregates, the filler fraction of the particles 
at the percolation point is much lower. Even in case 
the particles touch, the dc conductivity (CTdc) of these 
composites at high filling fractions turns out to be or- 
ders of magnitude lower than of the bulk material, as 
was recently illustrated for a particular crosslinked epoxy 
composite with filler particles of Phthalcon-11,'' Co ph- 
thalocyanine crystallites of 100 nm size, and explained 
by purely structural arguments.-^ 

When crystalline particles with a diameter of less than 
10 nm instead of 100 nm are used, the small size of the 
particles may impose another important restriction to 
the maximal possible composite conductivity, which is 
due to the density of states (DOS) involved in the dc 
conductivity through the network of particle contacts. 
Compared to larger crystallites, this DOS can be strongly 
reduced by the charging energy^^iii^i^i 

We show how those microscopic parameters, which 
govern the charge-transport process across many decades 
of length scales, can accurately and consistently be deter- 
mined by ac (alternating current) dielectric spectroscopy 
from a few Hz to infrared frequencies. In particular we 
can address the parameters for Mott variable-range hop- 



ping, for heterogeneity-induced enhanced ac response, 
for phonon- or photon-assisted nearest-neighbor hopping, 
and for the Drude response of individual nanocrystals. 
Due to these quantitative results we can unambiguously 
determine also the role of the nanocrystal charging en- 
ergy in limiting the hopping process. We apply the 
method to antimony-doped tin-oxide (ATO) crystallites 
of 7 nm diameter and to 100 nm sized crystallites of 
Phthalcon-11. It turns out that in densely packed crys- 
tallites of ATO, due to the strong infiuence of the charg- 
ing energy on the DOS, Udc at room temperature is four 
orders of magnitude lower than the dc conductivity ex- 
trapolated from the Drude plasma frequency (i^pu) of the 
crystallites - a result with obvious implications for the de- 
sign of conducting composites. The dielectric method is 
well suited for a variety of systems with restricted ge- 
ometries, as we will illustrate by a short discussion of 
phase-change materials^^ and granular oxidesJ^. 



II. CONDUCTIVITY IN CONDUCTING 
POLYMER COMPOSITES AND GRANULAR 
METALS 

For randomly placed conducting spheres in an insulat- 
ing matrix, the relation between adc and the fraction (jj 
of spheres is known from percolation t heory>i i ^ ^ i ^ ^ Above 
the percolation threshold a = ao\(l) — (l)c\* where (j)c « 0.16 
is the percolation threshold, t w 2.0 and ctq is approxi- 
mately equal to the conductivity of the spheres.™'^^ When 
the building blocks of the network are fractal aggregates 
instead of solid spheres, ctq has to be replaced by the 
aggregate conductivity (Ta and depends on the particle 
conductivity and, via the non-linear relation (Ta/cTp — 
(0p,c/O.16)^^^ with a; > 0, on the real percolation thresh- 
old 0p_c of the particles." The value of the exponent x is 
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related to the random- walk dimension and the fractal di- 
mension df, and is maximally 2(df — l)/(3 — df). This 
shows that on purely geometrical grounds for a network 
with (fip^c — 0.0055, at the highest filling fraction of the 
aggregates adc will be three to four orders of magnitude 
lower than in the pure filler powder.— 

As remarked in the introduction, when nanosized par- 
ticles are used as fillers, charging energies (and quantum 
size effects) impose a further important restriction to the 
maximal possible composite conductivityj^ii This effect 
can be conveniently studied in densely packed powders 
of filler material by dielectric spectroscopy. 

A. Parameters for the dc conductivity 

In the ohmic regime, if there is a non-negligible density 
of states around the chemical potential, and the temper- 
ature T is high enough that also the Coulomb interaction 
can be neglected (T > Tcrit), cr dc will obey Mott's equa- 
tion for conduction via variable-range hopping (VRH):^^ 

CTdc fx exp[-(ro,Mott/r)''] (I) 

with ly ~ 1/A and 

fcsTb.Mott ~ 20/(ghopa^), (2) 

where a denotes the decay length of the electron density, 
kB the Boltzmann constant, and ghop the density of states 
relevant in the hopping processJ^ 

For randomly packed spheres of radius R and spacing 
s the localization length a will be enlarged f^i^ i ^"^^^^ and 
can be approximated by^'^^ 

d = (2R/s)a. (3) 

In the following we drop the tilde. 

Below Tcrit the T dependence of the conductivity will 
be dominated by a soft Coulomb gap, leading to so-called 
Efros-Shklovskii (ES) VRH>ia 

adc cxexp[-(To,Es/r)'/']. (4) 

In the ES VRH model in the dilute limit of a large dis- 
tance between the particles Tq^es is given by 

To^ES = 2.8eV(47reeoaA:B), (5) 

with e the electron charge, eo the vacuum dielectric 
constant, and e the relative dielectric constant of the 
medium. Tent is given by T„it ~ e'*a5hop/[fcB(47reeo)^] 
and (for s > R) the charging energy by ^ e^/(47reeo_R). 

For densely packed small particles, at high temper- 
atures but still in the regime, where Coulomb interac- 
tions are important (T < Tcrit), ES VRH behavior will 
evolve into nearest-neighbor hopping at a temperature 
Ta. Above Ta the conduction is thermally activated with 
an activation energy AEa of the order of the charging 
energy, and Ta {AEAf/To^ES- 



The experiments of Yu et on thin films of highly 

monodispersed semiconducting nanocrystals of CdSe of 6 
nm diameter, slightly smaller than the ATO crystallites 
discussed here, showed good agreement between the the- 
oretical and experimental value of Tq^es, Tcrit and Ta-— 

B. Sub-THz and far-infrared regime 

At sufficiently low frequencies the conductivity will be 
frequency independent and equal to its dc value, because 
the inhomogeneities are averaged out by the motion of 
the charge carriers. The minimal length scale for homo- 
geneity is referred to as Thom, 

■^hom = (^dckBT/nhome^ifos/^D), (6) 

where nhom is the density of the carriers involved in the 
hopping process at the border of the homogenous regime, 
and the onset frequency fos for the frequency depen- 
dence of (J is divided by 2D, with D the dimension of 
the system. -'^^ 

At high enough frequencies, when during half a period 
of the oscillation of the applied field electrons can hop 
solely between nearest-neighbors, the major contribution 
to the conductivity will be due to tunnelling between 
localized states at neighboring sites (the pair limit) 
This incoherent process can be either by phonon-assisted 
or photon-assisted hopping, where in the latter case the 
energy difference between the sites is supplied by photons 
instead of phonons.^^ The phonon-assisted contribution 
to the conductivity is given by 

2 

crpho„o„(w) = ^e^wfcBTijj^pgg^pln^ {~^) ' ^''^ 

with Thop ^ a the decay length of the electronic state 
outside the conducting particles, ghop the relevant DOS 
at the Fermi energy Ep, and Wph the phonon 'attempt' 
frequencyii^ This formula is valid when to < Wph; at 
higher w, where the contribution of phonon-assisted hop- 
ping to a becomes constant, photon-assisted processes 
usually take over, with a conductivity CTphoton given by 

CTphoto„(t^) = ye^/tw^L^^pfffj^p In* (^j^^ ■ (8) 

The energy fceT in eq. ([7]) is in eq. ([8]) replaced by huj 
and the phonon attempt frequency ujph by 2/o, with Iq 
being the 'overlap' pre-factor for the energy levels of two 
neighboring sites. In analogy with t^ph, 

lo/h can be in- 
terpreted as the attempt frequency for photon-assisted 
hopping. Equation ([8|) is only valid when lo < lo/h. As 
in phonon-assisted hopping, a passes over into a plateau 
at high to. 

C. Visible and optical response 

At high frequencies (for ATO in the infrared regime) 
the short period of the electromagnetic field will restrict 
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the motion of the carriers to the nanocrystaUite, and the 
dielectric response characterized by the complex relative 
dielectric constant e{uj) = 1 — ujp-^/[uj{uj + iT)] will be 
Drude-like, with WpD the Drude plasma resonance fre- 
quency and r = l/r the damping rate. In practice the 
constant 1 has to be replaced by eoo due to other contri- 
butions in this frequency regime, like the polarization of 
the ion cores. ■^'^ The Drude plasma frequency is related 
to the number of carriers per unit of volume ncrys and 
the effective mass m* as 



-'pD 



2 / + 



(9) 



For damping rates comparable to the Drude plasma fre- 
quency, the real plasma frequency (where the dielectric 
constant becomes zero) will be larger than cOpD. T is 
determined by the boundaries of the nanoparticle and 
additional (ionized impurity) scattering: 



r = l/r = vf/L. 



crysi 



(10) 



where 1 /icrys is the sum of the inverse size of the particle 
and the inverse phonon scattering length. 
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FIG. 1: dc conductivity as function of T for a densely packed 
powder of 7%-doped ATO. For T < 50 K the data can be 
fitted by eq. (|4)| (a), while for T > 50 K the T dependence is 
activated (b). 



III. EXPERIMENTAL PROCEDURES AND 
DATA 

Measurements were performed on Sb-doped tin-oxide 
nanoparticles with [Sb]/([Sn]-|-[Sb]) equal to 0, 2, 5, 7, 
9, and 13 at.%. The particles are monocrystalline and 
spherical with diameters close to 7 nm.^"* The diameter 
of the 7% doped crystallite is 7.1 nm. Sb is incorporated 
in the casserite Sn02 lattice by replacing Sn^+. At the 
doping level of 7%, Sb is mainly present as Sb^"*", result- 
ing in n-type conductivity of the ATO particles according 
to Nutz et al.^ The amount of Sb'^+ present in the par- 
ticles is negligibly small i^i^ 

The followed experimental procedures for the dc con- 
ductivity and dielectric measurements are described in 
Refs. d and [l^. The thickness of the samples was typi- 
cally a few mm. The dc conductivity measurements were 
performed in the dark under helium atmosphere. 

The T dependence is given in Fig.[l] and the frequency 
dependence at temperatures down to 7 K in Fig. [51 All 
data shown are for 7%-doped ATO. Similar results were 
obtained at other doping levels, be it with different ab- 
solute values. The data were taken in the ohmic regime. 

The infrared (IR) transmittance was measured on a 
pellet of KBr mixed with a small amount of ATO. For 
the IR-reflectance we used a precipitated film of ATO 
with a thickness of about 1 mm. The data are shown 
in Fig. [21 For the analysis we also used the sub-THz 
transmittance and phase data (only shown in Fig. [?]). 

In Table [J we summarize our data on densely packed 7- 
nm-sized ATO crystallites and compare them with mea- 
surements on doped tin oxide published in the litera- 
ture. The values of agree within a factor 2, while 
the spread in the scattering rates is larger. 
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TABLE I: Transport parameters obtained for ATO and in- 
dium tin oxide (ITO) (second row) at doping levels of 10^" 
to lO^'^ per cm'^. The room-temperature dc conductivities are 
given in S/cm, the Drude frequencies and damping rate in 
s~^, and the effective mass in free-electron masses. The first 
and second row are obtained for films of AT02^ and ITO^^ 
resp., the last two rows contain the data on powders of 6% 
doped ATO particles of Nutz et al?^ and our data (labelled 
as pw for present work) on samples with 7% Sb doping. For a 
bulk material with a Drude frequency of 10^^ Hz, a scattering 
time of 10~^* s and a carrier mass of 0.3 mc a dc conductivity 
is expected of 10^ S/cm 



IV. ANALYSIS AND DISCUSSION 



In the analysis we first show the procedure to extract 
the parameter values from the data in the different fre- 
quency regimes and to check their consistency. We also 
make a comparison to the parameter values of Phthalcon- 
11, for which the data are published elsewhere Then 
we concentrate on the density of states; the latter being 
important for the dc conductivity. Subsequently, the im- 
plications for the use of the particles as fillers in nanocom- 
posites are discussed. 
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FIG. 2: Frequency dependence of the conductivity for 7%- 
doped ATO at various temperatures. In this double loga- 
rithmic plot the linear dependence at low temperatures is in 
agreement with eq. ((7}. 



1.0 



0.8 



0.6 



0.4 



0.2 



0.0 



Transmission (a) 



data 
■fit 



n ' 1 
,., Mil 



Reflectance (b) 

data 
- - fit 



100 200 300 400 100 200 300 400 

Energy (meV) 



FIG. 3: Transmittance and reflectance of ATO with 7% Sb 
doping in the infrared, after background correction, (a) For 
the transmission ATO is mixed with KBr and pressed into 
a 0.57 mm thick pellet to have sufficient transparency. The 
oscillations in the fit to the transmittance are an artefact (the 
fit is based on the calculated effective ATO-film thickness of 
0.005 mm), (b) Like for the sub-THz data, the reflection data 
are taken from precipitated films of 1 mm thick typically. The 
Drude fit is discussed in the text. The data compare well with 
the data published in the literature, see Table |T] 



A. Procedure 

Regarding the T dependence of a (Fig. [1]), the data 
can be fitted with v ~ 1 (eq. H])) if the fit is restricted to 
T > 50 K and with v = 0.5 (eq. (gl)) for T < 50 K. The 
I' = 1 fit gives an activation energy of 10^ K, while the 
exponent u = 0.5 at low T gives Tq^es = 0.6 • 10^ K. The 
localization length from Tq.es, see eq. ([S]), is calculated to 
be a = 3 nm. Using eq. ([3]) and (2i?/s) — 10 estimated 



from the packing density, we find a ^ 2 nm, in good 
agreement with the value calculated from Tq^es- 

The onset of the frequency dependence of the conduc- 
tivity (see Fig. [2]) signals that the carrier starts to feel 
the inhomogeneity of the underlying structure. Using 
eq. d?]), the typical length scale £hom at the onset can 
be found. For 7%-doped ATO at 300 K, the onset fre- 
quency /os = 3-10* Hz and CTdc — 10~^ S/cm give a value 

of ^hom"-hom = 9 • lO"^ Cm'T 

The linear frequency dependence of the conductivity 
at 7 K in the double logarithmic plot of Fig. [21 is in 
agreement with phonon-assisted tunneling, see eq. ([7|) . In 
the range of 10 - 100 cm"-'^ photon-assisted processes take 
pygj. 30 Applying eq. (O to the conductivity data at 293 
K and taking the usual value for the phonon frequency 
in solids = 10^2 s-ifi§.^ we find L^pffhop = 3 • 10^ 
eV^^cm^^, see Fig. |4l 

Turning next to the high-frequency data presented in 
Fig. [3l we performed a simple Drude analysis. The fit 
(tjpD — 11000 cm^^ and t = 3300 cm^^, together with 
a dielectric constant of 4.0) reproduces the main features 
of the increase of the transmission and the level of the 
refiectance (the oscillations in the fit to the transmittance 
are an artefact because the effective ATO film thickness 
of 0.005 mm is much smaller than the real thickness of the 
pressed KBr pellet). The number of carriers /idys of 10^^ 
cm~'^ is directly derived from the Drude frequency and is 
slightly lower than obtained from a simple interpretation 
of the chemical composition. The bulk dc conductivity 
calculated from the Drude parameters is 10^ S/cm. The 
fit parameters of the present samples are given in Table 
Uand agree well with the literature. 

Fig. |4] shows the reconstructed conductivity of ATO 
as function of frequency due to the processes discussed 
above. 

For ATO the important values for the dc conductivity 
can be deduced from the combination of variables that 
we found from the previous analysis (i) -^^hom^-hom = 9 • 

" -1 -r 3 • 10^ 



10^ 



(ii; 



i/hop^hop 



10^3 J-^m- 



eV ^cni ^ for photon-assisted hopping and lO''^ J ^m ^ 



''crys 



for the phonon-fit to the data at 7 K, and (iii) ric 
10^1 cm-3 and r = IQ-'^^ s. 

Using (in) the 'extrapolated' dc conductivity is 10^ 
S/cm, a factor 10* larger than the found value of 
10~2 S/cm. The estimated Fermi energy Ey is around 
2 eV, and from g{E) = (3/2)(n/£;F) (vahd for free elec- 
trons) we get giEp) — 5 • 10^^ eV^cm"^. From (ii) with 
^hop = a of 3 nm, we find for ghop = 3 • 10^^ eV~^cm~'^, 
a factor 10'^ lower than g{Ep). Note that this is an av- 
eraged density of states involved in photon assisted hop- 
ping. Due to the curvature of the density of states around 
the chemical potential, g{E) will be lower at lower ener- 
gies. For example, for the phonon-fit at 7 K g^^op is equal 
to 0.3 • 10^^ eV-^cm^^ 

The values for r and Lhom can be used as a consis- 
tency check. The combination of the estimated Fermi 
velocity of 0.8 • 10^ cm/s, with the crystallite size of 7 nm 
and m* = 0.3mc,^^ predicts a surface scattering rate of 
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FIG. 4: The various contributions to the conductivity as func- 
tion of photon energy. At room temperature at the lowest 
frequencies the conductivity is dominated by charging ener- 
gies, at intermediate frequencies phonon and photon assisted 
hopping processes describe the frequency dependence, and in 
the infrared the Drude conductivity inside the crystallites is 



are semiconducting crystals with a low number of charge 
carriers. All charges participating in the conductivity 
within the crystal also contribute to the dc conductiv- 
ity. As for ATO the obtained conduction parameters for 
Phthalcon-11 from the dielectric scans give a consistent 
picture. 



B. Density of states 

For ATO, the differences between the density of states 
involved in the hopping process ghop = 3 • 10^* eV^^cm""^ 
and the Drude conduction within the crystallites g{E-p) — 
5 • 10^^ eV~^cm~'^ are clearly significant. The result is 
as anticipated from the estimated charging energy of the 
order of 50 meV, and shows its importance for the dc 
powder conductivity. 

For Phthalcon-11 the very low number of carriers in- 
volved in the hopping process is similar to the number 
of carriers that determines the Drude contribution in the 
crystallites. Since the mean size of the particles is 20 
times larger than for ATO, the charging energies will be 
of the order of 3 meV, and hence are expected to be neg- 
ligible at room temperature. 



10"'^^ s~^, in agreement with the found value of F. Next, 
from ghopknT ~ nhom eV~^cm~^ we now can estimate 
'T-iiom at T = 293 K as rihom = 10^'' cm~'^. Using (i) and 
'T-hom we find I/hom = 0.3 /im. 

In short, the dielectric data of ATO allow a consis- 
tent picture of the conduction process. In these densely 
packed crystallites the localization length is enhanced 
by a factor 10 and the density of states involved in the 
dc conductivity is more than a factor 10'^ smaller than 
that in the conduction within the crystallites. The rel- 
atively large length scale for homogeneity is indicative 
for the presence of aggregates. Indeed, like in Ketj en- 
Black,^ nanoparticles of ATO are known to form chem- 
ically bonded aggregates that survive the preparation 
stagei^i Due to the nature of the chemical bond, the con- 
ductivities between neighboring crystallites in and out- 
side the aggregates are expected to be only slightly dif- 
ferent. Note that also the value of gays has to be seen 
as an average, as inhomogeneities in the doping of ATO 
might be present as welli^i^ 

For the studied Phthalcon-11 crystallites Wpo = 10^'^ 
s"-"^ and T — 10^^^ s leading to ricrys = 2 • 10^^ cm^^, 
i.e. about 1 charge per crystallite)^ The other values 
found for Phthalcon-11 are: (i) -Z^hom'^hom = 2-10^ cm~^, 
(ii) ffhop-^hop = 10^^ J"^m-i or 2.5 • 10^ eV-^cm-i. In 
these organic crystals with such a low carrier density, 
the charge carriers can be seen as an electron gas with 
an energy scale set by k^T, and g{E) can be estimated 
from g{E)kBT - ricrys to be 10^^ eV"i cm ^. This value 
of g{E) is the upper limit for t/hop and ghom- From ghom — 
10^^ eV~^cm~'^, we find a decay length a of 3 nm, as 
expected from the packing. 

The Phthalcon-11 parameters show that the crystals 



C. Implications 

In polymer nanocomposites with building blocks 
formed by diffusion-limited cluster aggregation, the airy 
structure of the particle network gives a strong reduc- 
tion in conductivity of the composite compared to the 
filler (for the Phthalconll /polymer composite a factor 
10^)»^ This effect can be compensated by using better 
conducting particles. Particles of ATO or ITO seem to 
be well-suited as the material is known to be very well- 
conducting. In addition, ATO crystallites are relatively 
easily obtained in sizes around 7 nm, and 'when properly 
dispersed can give polymer composites with a low per- 
colation thresholdi^ However, even if the filler nanopar- 
ticles in the composite touch, they will not be in better 
contact than in a densely packed powder. As shown here 
for ATO, for these small crystallites the DOS involved in 
(Tdc is dramatically reduced due to the shift of the energy 
levels away from the Fermi level by Coulomb charging 
effects. As a consequence, an additional four orders of 
magnitude in a^c are lost compared to the bulk value. 

Other systems where size restrictions are expected to 
be present might be conveniently studied in a similar 
way. For example several chalcogenide alloys exhibit 
a pronounced contrast between the optical absorption 
in the metastable rocksalt after the intense laser pulse 
and in the initial amorphous phase. As shown by ex- 
tended x-ray absorption fine structure spectroscopy (EX- 
AFS) the resistive change after the intense laser recording 
pulse goes together with a crystallization process, where 
also small domains are inherently present. Our dielectric 
method might visualize to what extent the domain walls 
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after crystallization limit the conductivity and have con- 
sequences for the band structure calculations. If the walls 
become real barriers quantum size effects in the small do- 
mains will invalidate the use of periodic boundary condi- 
tions in the calculations. Also the glassy behavior in the 
conductance of deposited indium-tin oxide samples in the 
insulating regime)^ and of quench-condensed insulating 
granular metals'^^ might be further clarified by the use of 
our dielectric approach and analysis. Scanning the fre- 
quency will reveal the evolution of the length scales and 
DOS involved in the relaxation processes. 

V. CONCLUSIONS 

By combining data of sub-THz transmission with in- 
frared transmission and reflection we were able to ex- 
plain the full frequency response of densely-packed nano- 
sized crystallites using the parameters for Mott variable- 
range hopping, for heterogeneity-induced enhanced ac re- 



sponse, for phonon- or photon-assisted nearest-neighbor 
hopping, and for the Drude response of individual 
nanocrystals. For 7 nm antimony-doped tin-oxide parti- 
cles the analysis unambiguously quantified the reduction 
of the density of states involved in the dc conduction 
compared to the value extrapolated from the Drude re- 
sponse at infrared frequencies. Dielectric scans with a 
similar analysis will also be revealing in other systems 
where size limitations are expected to play a role. 
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